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The  ability  of  two-equation  turbulence  models  to  accurately  predict  separated  flows  is  ana¬ 
lyzed  from  a  combined  theoretical  and  computational  standpoint.  Turbulent  flow  past  a  backward 
facing  step  is  chosen  as  a  test  case  in  an  effort  to  resolve  the  variety  of  conflicting  results  that  have 
been  published  during  the  past  decade  concerning  the  performance  of  two-equation  models.  It  is 
found  that  the  errors  in  the  reported  predictions  of  the  K- E  model  have  two  major  origins:  (1)  nu¬ 
merical  problems  arising  from  inadequate  resolution,  and  (2)  inaccurate  predictions  for  normal 
Reynolds  stress  differences  arising  from  the  use  of  an  isotropic  eddy  viscosity.  Inadequacies  in 
near  wall  modeling  play  a  substantially  smaller  role.  Detailed  calculations  are  presented  which 
strongly  indicate  that  the  standard  K-t  model  —  when  modified  with  an  independently  calibrated 
anisotropic  eddy  viscosity  —  can  yield  surprisingly  good  predictions  for  the  backstep  problem. 


+  Permanent  Address:  Department  of  Mechanical  Engineering,  Stevens  Institute  of  Technology,  Hobo¬ 
ken,  NJ  07030 

*  This  research  was  supported  by  the  National  Aeronautics  and  Space  Administration  under  NASA  Con¬ 
tract  No.  NAS1-18605  while  the  authors  were  in  residence  at  the  Institute  for  Computer  Applications  in 
Science  and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23665. 


1.  INTRODUCTION 


Turbulent  flow  past  a  backward-facing  step  has  played  a  central  role  in  benchmarking  the  per¬ 
formance  of  turbulence  models  for  separated  flows.  During  the  past  decade  —  beginning  with  the 

1980/81  Stanford  conference  on  complex  turbulent  flows1  —  a  variety  of  two-equation  turbulence 
models  were  tested  and  compared  with  the  experimental  data  of  Kim,  Kline  and  Johnston^  and 

Eaton  and  Johnston3  for  the  backstep  problem.  Initial  results 1  indicated  that  the  standard  K-z  mod¬ 
el,  with  wall  functions,  underpredicted  the  reattachment  point  by  a  substantial  amount  of  the  order 
of  20-25%.  A  significant  number  of  studies  have  been  subsequently  published  using  alternative 
forms  of  the  K-z  model  wherein  a  variety  of  conflicting  results  have  been  reported.  For  example, 

Sindir4  made  modifications  to  account  for  streamline  curvature  based  on  the  algebraic  stress  model 

of  Gibson5;  some  improvements  were  obtained  although  the  results  were  somewhat  mixed.  Chen6 
performed  calculations  with  a  multiple  scale  K-z  model  wherein  the  turbulent  kinetic  energy  K  and 
the  turbulent  dissipation  rate  e  were  decomposed  into  low  and  high  wavenumber  parts  along  the 

lines  suggested  by  Hanjalic  et  al?  Significantly  improved  results  for  the  reattachment  point  were 
obtained  (the  underprediction  was  reduced  to  5%),  however,  no  detailed  comparison  of  the  Rey¬ 
nolds  stresses  were  made.  Speziale8  reported  comparable  improvements  for  the  backstep  problem 
based  on  an  anisotropic  K-z  model  which  appeared  to  suggest  that  the  main  source  of  the  errors 

could  be  due  to  the  use  of  an  isotropic  eddy-viscosity  in  the  standard  K-z  model.  Later,  So  et  al 9 
presented  results  which  seemed  to  indicate  that  the  right  near  wall  model  could  significantly  im¬ 
prove  the  predictions  of  the  K-z  model  for  the  backstep  problem.  Similar  claims  were  made  by 

Kamiadakis  et  a/10  who  argued  that  the  new  near  wall  treatment  in  their  RNG  based  K-z  model 
gave  rise  to  drastically  improved  predictions  for  the  reattachment  point.  However,  Avva,  Kline 

and  Ferziger1 1  presented  results  which  suggested  that  the  large  underprediction  of  the  reattachment 
point  attributed  to  the  standard  K-z  model  was  due  mainly  to  the  fact  that  the  previously  reported 
computations  were  under-resolved.  When  sufficient  resolution  was  used,  they  found  that  the  actual 
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error  was  only  of  the  order  of  10%  —  a  deficiency  which  they  eliminated  by  introducing  an  empir¬ 
ical  correlation  for  streamline  curvature.  Considering  the  need  to  accurately  predict  separated  tur¬ 
bulent  flows  —  which  can  have  a  wealth  of  important  scientific  and  engineering  applications  —  it 
is  rather  unsetding  that  such  a  wide  range  of  conflicting  claims  still  permeate  the  literature.  This 
establishes  the  motivation  of  the  present  paper  which  is  to  attempt  to  clarify  this  issue. 

In  this  paper,  the  following  questions  will  be  addressed: 

1)  Precisely  what  does  the  standard  K-e  model  predict  for  this  backstep  problem  when  the  numer¬ 
ics  are  done  properly  to  insure  full  resolution  of  the  flowfield  and  accurate  implementation  of 
the  boundary  conditions, 

2)  What  is  the  physical  source  of  actual  errors  in  the  standard  K-e  model,  and 

3)  Can  these  errors  be  eliminated  within  the  framework  of  two-equation  turbulence  models  with¬ 
out  the  ad  hoc  adjustment  of  constants  or  the  introduction  of  other  arbitrary  empiricisms? 

These  questions  will  be  addressed  by  conducting  sufficiently  resolved  computations  of  the  stan¬ 
dard  and  anisotropic  K-e  model  —  with  several  different  wall  boundary  conditions  —  for  the  back- 

step  configuration  of  Kim  et  al ?  and  Eaton  and  Johnston3.  By  making  detailed  comparisons  with 
the  experiments,  it  will  be  demonstrated  that  the  standard  K-e  model  yields  results  that  are  within 
12%  of  the  experimental  data.  With  the  introduction  of  an  anisotropic  eddy  viscosity,  calibrated 
independently  without  any  ad  hoc  empiricisms,  the  K-e  model  is  shown  to  yield  an  exceptionally 
good  prediction  for  the  reattachment  point  that  is  within  a  few  percent  of  the  experimental  result. 
The  physical  implications  that  these  results  have  on  the  previous  work  cited  above,  as  well  as  on 
future  applications  to  more  complex  separated  flows,  will  be  discussed  in  detail  in  the  sections  to 
follow. 
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2.  FORMULATION  OF  THE  PHYSICAL  PROBLEM 


The  problem  to  be  considered  is  the  fully-developed  turbulent  flow  of  an  incompressible  vis¬ 
cous  fluid  past  a  backward-facing  step  (a  schematic  is  provided  in  figure  1).  Calculations  will  be 

conducted  for  the  Kim  et  al 2  configuration  wherein  the  expansion  ratio  (step  height:  outlet  channel 
height)  £  is  1:3  and  the  Reynolds  number  Re  =  132,000  based  on  the  inlet  centerline  mean  velocity 
and  outlet  channel  height.  The  Reynolds  averaged  Navier-Stokes  and  continuity  equations  are 
solved  which  take  the  form: 
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where,  u  and  v  are  the  mean  velocity  components  in  the  x  andy  directions;  p  is  the  modified  mean 


pressure;  t^,  x^  and  x  are  the  components  of  the  Reynolds  stress  tensor  x-  s  u\u'j\  and  v  is  the 
kinematic  viscosity.  In  the  standard  K-e  model  with  isotropic  eddy-viscosity  the  Reynolds  stress 
tensor  takes  the  form  (see  Launder  and  Spalding12) 
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is  the  mean  rate  of  strain  tensor,  K=  -xu  is  the  turbulent  kinetic  energy,  e  is  the  turbulent  dissi¬ 
pation  rate,  ui  =  (u,  v)  is  the  mean  velocity  vector,  and  is  a  dimensionless  constant  which  is 
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taken  to  be  0.09.  An  anisotropic  K-e  model  will  be  considered — namely,  the  nonlinear  K-e  mod¬ 
el  of  Speziale8  —  wherein  the  Reynolds  stress  tensor  takes  the  form 


\j  =  |K8..  -  2C^Stj  -  4CdC^  <4  -  +  SikSkJ  -|sk,Sw 8y)  (6) 


where 
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dt  *dxk  Bxk  *J  dxk 
is  the  frame-indifferent  Oldroyd  derivative  and  Cq  is  a  dimensionless  constant  which  takes  on  the 
value  of  1.68  based  on  a  calibration  with  turbulent  channel  flow  data.  The  standard  K-e  model  is 
recovered  in  the  limit  as  CD  — >  0.  Both  the  standard  and  anisotropic  K-e  model  are  solved  in  con¬ 
junction  with  modeled  transport  equations  for  the  turbulent  kinetic  energy  K  and  turbulent  dissipa¬ 
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where  vr  s  —  is  the  eddy  viscosity, 
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is  the  turbulence  production,  and  Cel,  C&,  and  oe  are  dimensionless  constants  which,  respec¬ 
tively,  are  taken  to  be  1.44, 1.92, 1.0  and  1.3.  The  Reynolds  averaged  equations  (l)-(7)  are  solved 
subject  to  the  following  boundary  conditions: 

(a)  inlet  profiles  for  u,  K  and  e  are  specified  five  step  heights  upstream  of  the  step  comer  (ti  is  taken 
from  the  experimental  data3  and  the  corresponding  profiles  for  K  and  e  are  computed  from  the 
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model  formulated  for  channel  flow), 

(b)  the  law  of  the  wall  is  used  at  the  upper  and  the  lower  walls,  and 

(c)  standard  extrapolated  outflow  conditions  are  applied  thirty  step  heights  downstream  of  the  step 
comer. 

The  law  of  the  wall  is  applied  in  two  different  forms.  In  the  standard  two-layer  form  of  the  law 
of  the  wall: 

-+  1  ,  +  .  K  -,-1/2  ~3/4tf3/2  .... 

u  =  -  Iny  +5,  —  =  ,  e  =  — —  (10) 

K  ux  y 

at  the  first  grid  pointy  away  from  the  wall  if  y  +  =yux/v  >  11.6  given  that  u  +  =  u/ux  (i ^  is  the 

shear  velocity  and  k  =  0.41  is  the  von  K£rm£n  constant);  if  y  +  <  11.6,  then  u,  K  and  e  are  inter¬ 
polated  to  their  wall  values  based  on  viscous  sublayer  constraints.13  A  three-layer  law  of  the  wall 
where 

y+,  fory+£5 

-3.05+5  lny+,  for5<y+£30  t11) 

5.5  +  2.5  lny+,  fory+>30 

is  also  implemented  along  the  lines  of  Avva  et  al. 1 1  wherein  the  fact  that  the  normal  derivative  of 
K  vanishes  at  the  wall  is  made  use  of  along  with  more  elaborate  interpolation  formulas  for  K  and 

e  (see,  also,  Amano13  for  a  detailed  discussion).  It  must  be  said  at  the  outset  that  the  law  of  the 
wall  does  not  formally  apply  to  separated  turbulent  boundary  layers.  However,  since  the  separation 
point  is  fixed  at  the  comer  of  the  backstep  —  and  the  flowfield  is  solved  iteratively  so  that  the  shear 
velocity  ttj  can  be  updated  until  it  converges  —  major  errors  do  not  appear  to  result  from  its  use  as 
will  be  demonstrated  in  the  next  section. 

The  governing  equations  (l)-(7)  are  solved  using  a  finite  volume  method  (see  Lilley  and 
Rhode14  and  Thangam  and  Hur15).  Here  we  are  interested  in  the  steady  state  solution  which  is  ob- 
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tained  by  solving  the  system  of  algebraic  equations  by  a  line  relaxation  method14  with  the  repeated 
application  of  the  tridiagonal  matrix  solution  algorithm.16  The  issue  of  resolution  is  crucial  for  the 
backstep  problem.  As  discussed  by  Avva  et  a/.,11  several  researchers  have  reported  results  that 

were  in  severe  error  due  to  a  lack  of  adequate  resolution.  Thangam  and  Hur15  conducted  a  careful 
grid  refinement  study  based  on  this  finite  volume  method  for  grids  containing  150x75  to  400x200 
mesh  points  (these  meshes  have  variable  grid  spacing  to  allow  for  maximum  resolution  near  the 
step  comer  and  the  walls).  The  conclusion  of  their  study  was  that  a  200x100  mesh  yielded  results 
that  were  within  0.3%  of  the  grid  independent  solution;  the  use  of  significantly  coarser  meshes 
could  lead  to  appreciable  errors  whereas,  for  the  most  part,  the  improvement  to  be  gained  by  the 
use  of  finer  meshes  would  be  minimal.  All  of  the  computations  conducted  in  this  study  were  per¬ 
formed  using  this  200x100  nonuniform  mesh.  As  indicated  earlier,  the  inlet  conditions  were  spec¬ 
ified  5  step  heights  upstream  of  the  step  corner  and  the  outlet  boundary  conditions  were  specified 
30  step  heights  downstream  of  the  step  comer.  It  is  crucial  that  a  sufficient  distance  downstream 
of  the  reattachment  point  be  allowed  before  the  outflow  conditions  are  imposed.  Many  earlier 
computations  of  the  backstep  problem  were  in  significant  error  due  to  the  imposition  of  fully-de¬ 
veloped  outflow  conditions  too  close  to  the  reattachment  point. 

The  steady  state  solution  of  Eqs.  (l)-(7)  is  obtained  by  an  iterative  solution  of  the  discretized 
equations.  The  computed  solution  was  assumed  to  have  converged  to  its  steady  state  when  the  root 
mean  square  of  the  average  difference  between  successive  iterations  was  less  than  10"4  for  the  mass 

source.14,15  Approximately  2000  iterations  were  needed  for  the  convergence  of  the  standard  K-z 
model;  this  corresponds  to  approximately  30  minutes  of  CPU  time  in  a  partially  vectorized  mode 
on  the  CRAY-2S  supercomputer  using  64-bit  precision.  The  anisotropic  K-z  model  requires  ap¬ 
proximately  30%  more  CPU  time  due  to  the  fact  that  the  additional  correction  terms  to  the  standard 
K-z  model  have  to  be  evaluated  during  each  iteration.  These  correction  terms  are  dispersive  rather 
than  being  dissipative  —  an  additional  feature  that  slows  convergence. 
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3.  DISCUSSION  OF  THE  RESULTS 

First,  results  will  be  presented  for  the  standard  K-e  model  using  the  two-layer  law  of  the  wall 
boundary  conditions.  For  this  case  —  as  well  as  the  other  results  to  follow  —  computed  results  for 
the  mean  velocity  streamlines,  the  streamwise  mean  velocity  profiles,  the  streamwise  turbulence 

intensity  profiles  and  the  turbulence  shear  stress  profiles  will  be  compared  with  the  Kim  et  al ?  ex- 

'y 

perimental  data  as  updated  by  Eaton  and  Johnston.  In  figure  2(a)  the  computed  streamlines  are 
shown  which  indicate  reattachment  at  XIH  =  6.0  —  a  result  which  is  approximately  a  15%  under¬ 
prediction  of  the  experimental  reattachment  point  of  XIH  =  7.1.  This  is  in  contrast  to  earlier  report¬ 
ed  results  which  seemed  to  indicate  that  the  standard  K-e  model  underpredicts  the  reattachment 
point  by  20-25%  —  an  exaggerated  underprediction  due  to  the  lack  of  adequate  resolution  in  those 
computations.  (A  spurious  underprediction  of  the  reattachment  point  can  also  result  from  the  ap¬ 
plication  of  outflow  boundary  conditions  at  a  distance  XIH  <  25  from  the  step  comer  as  some  au¬ 
thors  have  done).  In  figure  2(b),  the  streamwise  mean  velocity  profiles  predicted  by  the  standard 
K-e  model  are  compared  with  the  experimental  data.  Except  in  the  vicinity  of  the  reattachment 
point,  the  comparisons  are  fairly  good.  More  serious  discrepancies  between  the  model  predictions 
and  the  experimental  data  occur  in  the  initial  part  of  the  recovery  zone  for  the  streamwise  turbu¬ 
lence  intensity  profiles  as  shown  in  figure  3(a).  However,  the  model  predictions  for  the  turbulence 
shear  stress  profiles  are  reasonably  good  as  can  be  seen  from  figure  3(b). 

Next,  we  will  examine  the  effect  of  the  near  wall  modeling.  In  figures  4(a)-(b),  the  mean  flow 
streamlines  and  streamwise  mean  velocity  profiles  obtained  from  the  standard  K-e  model  with  the 
three  layer  law  of  the  wall  boundary  condition  is  shown.  With  this  more  accurate  version  of  the 
law  of  the  wall,  the  standard  K-e  model  predicts  reattachment  at  XIH  =  6.25  —  a  result  which  is 

extremely  close  to  that  reported  by  Avva  et  al} 1  By  integrating  the  K-e  model  directly  to  the  wall 
—  with  an  asymptotically  consistent  low  Reynolds  number  correction  (cf.,  So  et  al .9  and  Speziale 

et  al. 17)  —  it  is  possible  to  improve  the  prediction  for  the  reattachment  point  by  approximately  5%. 
Some  earlier  published  work  claiming  that  the  near  wall  modeling  could  make  a  larger  difference 
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were  in  all  likelihood  due  to  inadequate  resolution  of  the  computations  (e.g.,  So  et  al 9  used  a 
93x66  mesh).  Considering  the  current  state  of  development  of  near  wall  modeling  —  as  well  as 
the  relatively  small  benefits  that  it  can  yield  for  the  backstep  problem  —  the  remaining  results  to 
be  shown  will  be  based  on  the  three-layer  law  of  the  wall  boundary  condition.  However,  for  tur¬ 
bulent  flows  where  the  separation  point  is  not  fixed,  near  wall  turbulence  modeling  can  play  a  much 
more  important  role. 

Now,  it  will  be  demonstrated  that  the  use  of  the  anisotropic  eddy-viscosity  model  (6)  can  yield 
a  more  significant  improvement  in  the  results.  The  computed  mean  velocity  streamlines  obtained 
from  the  nonlinear  K-e  model  are  shown  in  figure  6(a).  They  indicate  reattachment  at  X/H  ~  6.9: 
a  result  that  is  within  3%  of  the  mean  experimental  value  of  X/H  =  7.1.  This  more  accurate  predic¬ 
tion  of  the  reattachment  point  is  reflected  in  the  streamwise  mean  velocity  profiles  shown  in  figure 

6(b)  which  are  in  better  agreement  with  the  experimental  data2,3  than  those  obtained  from  the  stan¬ 
dard  K-e  model  based  on  an  isotropic  eddy  viscosity.  Some  of  the  earlier  reported  results  for  the 

nonlinear  K-e  model  (see  Speziale8  and  Speziale  and  Ngo18)  were  lower  due  to  insufficient  numer¬ 
ical  resolution  of  the  flowfield.  The  corresponding  streamwise  turbulence  intensity  profiles  and 
turbulence  shear  stress  profiles  are  compared  with  the  experimental  data  in  figure  7(a)-(b).  On  bal¬ 
ance,  the  agreement  between  the  model  predictions  and  experimental  data  is  good  (it  must  be  re¬ 
membered  that  the  uncertainty  in  the  turbulence  Reynolds  stress  measurements  is  of  the  order  of 
10%).  The  most  notable  difference  between  the  predictions  of  the  anisotropic  and  the  standard  K- 
e  model  lies  in  the  streamwise  turbulence  intensity  located  approximately  one  step  height  above 
the  comer  of  the  step.  Unlike  the  standard  K-e  model,  the  nonlinear  K-e  model  predicts  a  dramatic 
trough  in  this  region  —  a  result  that  is  consistent  with  the  more  recent  independent  experi¬ 
ments.19,20  In  fact,  it  is  the  more  accurate  prediction  of  the  normal  Reynolds  stress  difference 
x  - 1  that  accounts  for  the  better  predictions  of  the  nonlinear  K-e  model.  This  can  be  illustrated 

yy  xx 

directly  in  turbulent  channel  flow.  In  figure  8,  the  normal  Reynolds  stress  difference  predicted  by 
the  nonlinear  K-e  model  is  compared  with  the  experimental  data  of  Laufer.  It  is  clear  that  the 
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nonlinear  K-z  model  does  a  reasonably  good  job  of  reproducing  the  experimental  data  whereas  the 
standard  K-z  model  erroneously  predicts  that  %  —  %xx  =  0.  The  accurate  prediction  of  this  nor¬ 
mal  Reynolds  stress  difference  is  important  in  the  backstep  problem  since  the  mean  flow  stream- 
function  \j7  (where  u  =  -d\\i/dy  ,  and  v  s  d\\f/dx)  is  a  solution  of  the  transport  equation18 


D  d 

u— (V2\j/) +v^  (V2y)  =  vV4y- 


02(x  -x  )  d2x  02x 

'  4*4/  VV  -  V4I 


yy  xx 


xy 


xy 


dxdy 


dx 2 


dy4 


(12) 


and,  hence,  the  normal  Reynolds  stress  difference  x  -  x  contributes  directly  to  the  calculation 

yy  xx 

of  the  mean  velocity  field.  This  is  the  case  for  any  flow  with  streamline  curvature  since  the  normal 
Reynolds  stress  term  on  the  r.h.s.  of  (12)  can  be  converted  to  the  alternate  form 


d2 (x  -  x  ) 

'  MM  OO ' 


nn  ss ' 


dnds 


in  terms  of  an  intrinsic  coordinate  system  where  n  and  s  are,  respectively,  the  normal  and  tangential 
directions  along  a  streamline. 

Finally,  we  will  examine  the  RNG  based  K-z  model  of  Yakhot  and  Orszag.22  For  high  turbu¬ 
lence  Reynolds  numbers,  it  is  of  the  same  general  form  as  the  standard  K-z  model  given  by  Eqs. 
(4),  (8)  and  (9).  However,  the  numerical  values  of  the  constants  —  which  are  computed  by  the 
RNG  approach  —  take  on  the  alternate  values: 

C^=  0.0837,  Cel  =  1.063,  Ce2=  1.7215 
Og  =  0.7179,  oE  =  0.7179,  K  =  0.372 

For  the  low  turbulence  Reynolds  numbers  which  occur  close  to  the  walls  in  this  backstep  problem, 
the  RNG  K-z  model  has  a  built-in  correction  that  allows  for  a  direct  integration  to  a  solid  boundary 
with  the  no-slip  condition  applied.  However,  we  found  this  near  wall  correction  to  be  rather  am¬ 
biguous  —  especially  in  the  presence  of  geometrical  discontinuities  such  as  those  which  occur  near 
the  backstep  comers.  Hence,  we  will  present  results  in  which  we  match  to  the  three-layer  law  of 
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the  wall.  In  figures  9(a)-(b),  the  mean  flow  streamlines  and  the  streamwise  mean  velocity  profiles 
obtained  from  the  RNG  K-e  model  are  shown.  The  RNG  K-e  model  predicts  reattachment  at 
X/H  ~  4  —  a  result  that  constitutes  a  substantial  underprediction  of  the  experimental  reattachment 
point  of  X/H  =  7.1.  Due  to  this  significant  underprediction  of  the  reattachment  point,  the  stream- 
wise  turbulence  intensity  and  turbulence  shear  stress  profiles  predicted  by  the  RNG  K-e  model  are 
in  serious  error  as  illustrated  in  figures  10(a)-(b).  The  origin  of  this  substantial  underprediction  of 
the  reattachment  point  lies  in  the  choice  of  Cel  -  l  .063.  In  a  homogeneous  shear  flow,  it  can  be 
shown  that  the  eddy  viscosity 

\T~-exp{Xt*)  (13) 


where  t*  is  the  time  (nondimensionalized  by  the  shear  rate)  and  X  is  the  growth  rate  give  jy23 

,  2  -i  1/2 


X  = 


Cu(Ce2-  CEl) 


£lJ 


(Cel-1)  (Ce2-1) 


(14) 


Hence  the  growth  rate  of  the  eddy  viscosity  becomes  singular  as  Ce\  -» 1.  For  CEl  =  1.063  —  as 
opposed  to  the  more  standard  value  of  Cel  =  1.44  —  the  growth  rate  of  the  eddy  viscosity  will  be 
overly  large.  An  overprediction  of  the  eddy  viscosity  will  make  the  model  too  dissipative  —  a  fea¬ 
ture  which  will  cause  the  separation  bubble  to  shorten  as  shown  in  figure  9(a).  Hence,  while  the 
use  of  an  alternate  near  wall  treatment  with  the  RNG  K-e  model  could  cause  the  separation  bubble 
to  become  slightly  larger,  we  feel  that  it  is  very  unlikely  that  a  reattachment  point  close  to  X/H  ~  7 

could  be  obtained  from  this  RNG  K-e  model  as  reported  in  Kamiadakis  et  a/.10  (We  believe  that 
the  calculations  of  Karniadakis  et  a/.10  may  have  been  under- resolved).  It  should  be  noted,  how¬ 
ever,  that  Smith  and  Yakhot24  have  recently  determined  that  there  may  have  been  an  error  in  the 
original  calculation  of  Cei  by  Yakhot  and  Orszag22;  they  obtained  an  alternate  value  of  Ct\  =  1.40. 
With  this  higher  value  of  Cei  —  and  with  the  use  of  an  RNG  based  anisotropic  eddy-viscosity  (see 

Rubinstein  and  Barton25)  —  there  is  no  question  that  the  RNG  K-e  model  would  yield  excellent 
results  for  the  backstep  problem. 
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4.  CONCLUSIONS 


A  reasonably  detailed  theoretical  and  computational  analysis  of  the  performance  of  two-equa¬ 
tion  turbulence  models  for  turbulent  flow  past  a  backward  facing  step  has  been  presented  in  an  ef¬ 
fort  to  resolve  numerous  conflicting  studies  that  have  been  published  during  the  past  ten  years.  The 
following  general  conclusions  were  arrived  at: 

1)  The  standard  K-e  model  with  three-layer  wall  functions  predicts  reattachment  at  XI H  ~  6.25  — 
a  result  that  is  within  12%  of  the  experimental  value  of  Kim,  Kline  and  Johnston  for  this  back- 

step  problem.  This  result  is  in  close  agreement  with  the  recent  computations  of  Avva  et  a/.11 
and  hence  there  is  little  doubt  that  earlier  studies  which  attributed  a  20-25%  underprediction  of 
the  reattachment  point  to  the  standard  K-e  model  were  in  error  due  to  inadequate  numerical  res¬ 
olution  of  the  flowfield. 

2)  When  the  standard  K-e  model  is  appropriately  modified  to  include  an  anisotropic  eddy  viscos¬ 
ity,8  the  reattachment  point  moves  out  to  X/H  »  6.9.  This  result,  as  well  as  the  other  mean  ve¬ 
locity  statistics,  are  within  3%  of  the  experimental  data  —  an  improvement  due  to  the  more  ac¬ 
curate  prediction  of  normal  Reynolds  stress  differences. 

3)  The  RNG  K-e  model  of  Yakhot  and  Orszag22  substantially  underpredicts  the  reattachment 
point.  This  problem  arises  since  Cej  is  too  close  to  1  —  a  feature  that  causes  the  model  to  be 
overly  dissipative  (e.g.,  when  C£\  -  1,  the  growth  rate  of  the  eddy  viscosity  becomes  singular 

in  homogeneous  shear  flow).  The  previously  published  results  of  Kamiadakis  et  a/.,10  to  the 
contrary,  appear  to  have  been  under-resolved.  However,  there  is  no  question  that  with  the  re¬ 
cent  suggested  correction  (Smith  and  Yakhot24)  of  Cel  to  1 .40  —  and  with  the  use  of  an  aniso¬ 
tropic  eddy  viscosity  —  the  RNG  K-e  model  would  yield  excellent  results  for  the  backstep 
problem. 

Finally,  some  remarks  are  in  order  concerning  the  implications  that  these  results  have  on  tur¬ 
bulence  modeling.  The  deficiencies  of  two-equation  models  are  well  established,  particularly  in 

turbulent  flows  with  body  forces  or  Reynolds  stress  relaxation  effects.2^  Consequently,  the  find- 
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ings  of  this  study  should  not  be  interpreted  as  an  unequivocal  endorsement  of  two-equation  models. 
Nonetheless,  this  study  certainly  does  indicate  that  properly  calibrated  two-equation  turbulence 
models,  with  an  anisotropic  eddy  viscosity,  can  yield  excellent  results  for  the  backstep  problem  that 
are  far  superior  to  the  results  obtained  from  the  older  zero  or  one  equation  models. 
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(b)  Dimensionless  Mean  Velocity  Profile 
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Figure  6 
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Comparison  of  the  computed  values  of  the  normal  Reynolds  stress  dif¬ 
ference  obtained  from  the  nonlinear  K-e  model  with  the  experimental 
data  of  Laufer  (1951)  for  the  fully-developed  turbulent  channel  flow.  (It 
should  be  noted  that  the  standard  K- e  model  predicts  a  zero  normal 
stress  difference) 

Figure  8 
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layer  model;  o  Experiments  of  Kim  etal,  1980;  Eaton  &  Johnston,  1981] 

Figure  10 
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